Induced local spin-singlet amplitude and pseudogap in high Tc cuprates 
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In this paper we show that local spin-singlet amplitude with d-wave symmetry, (|Ad(0)p), can be 
induced by short-range spin correlations even in the absence of pairing interactions. Fluctuation the- 
ory is formulated to make connection between pseudogap temperature T* , pseudogap size Apg and 
(|Ad(0)|^). In the present scenario for the pseudogap, the normal state pseudogap is caused by the 
induced local spin-singiet amplitude due to short-range spin correlations, which compete in the low 
energy sector with superconducting correlations to make Tc go to zero near half-filling. Calculated 
T* falls from a high value onto the Tc line and closely follows mean- field Neel temperature T^^^ . 
The calculated Apg is in good agreement with experimental results. We propose an experiment in 
which the present scenario can be critically tested. 

PACS numbers: 71.10.Fd, 71.27.-|-a 



The recent discovery of pseudogap in underdoped high 
Tc cuprates has challenged condensed matter physicists 
for several years. The pseudogap behavior 0] is observed 
as strong suppression of low frequency spectral weight 
below some characteristic temperature T* higher than 
transition temperature Tc- This anomalous phenomenon 
has been observed in angle resolved photoemission spec- 
troscopy (ARPES), ||,| specific heat, g tmmeling, |] 
NMR, ||] and optical conductivity. |Q One of the most 
puzzling questions in pseudogap phenomena is why T* 
has a completely different doping dependence from Tc, in 
spite of their possibly close relation. 

Among several theoretical proposals |l| to understand 
the pseudogap phenomena, the superconducting (SC) 
fluctuation scenario |^ of pseudogap has received much 
attention from the physics community. This is because 
some experiments such as ARPES and tunneling exper- 
iments show that the normal state pseudogap has the 
same angular dependence and magnitude as the SC gap 
and that often the only difference between the spec- 
tra in the pseudogap state and the SC state is in their 
linewidths. The basic idea of this scenario is that SC gap 
amplitude forms at T^^^ while its phase coherence is es- 
tablished at Tc lower than T^^^ . Hence T* ~ T^^ and 
below this temperature SC fluctuations become stronger 
until they diverge at Tc. 

In spite of its success in explaining some features of the 
pseudogap, it suffers from at least three important draw- 
backs which were often overlooked in the past. First, 
just below T* there is no experimental evidence of char- 
acteristic features associated with SC fluctuations such 
as fluctuating diamagnetic (Meissner) effect, fluctuating 
superfluid density and so on. It appears that T* has noth- 
ing to do with superfluid "rigidity" . One experiment to 
strongly support this argument was recently carried out 
by Orenstein's group |9|. In their high-frequency conduc- 



tivity measurements tracking the phase-coherence time r 
in the normal state, the temperature Tq where the phase- 
stiffness of superfluidity disappears is at most 25 K above 
Tc for underdoped cuprates. Second, when SC correla- 
tions are treated on equal footing with antiferromagnetic 
(AF) correlations, which is more realistic from both theo- 
retical and experimental points of view, T^^^ goes down 
to zero with decreasing doping due to the competition 
with the AF correlations, as shown in Fig. |l|(a). Then 
the above scenario (T* ~ T^^^) is inconsistent even with 
the doping dependence of T*, which increases with de- 
creasing doping. Apparently experimentally observed T* 
stays in between T^^^ and T^^^ . Furthermore, in this 
situation the origin of the pseudogap itself is question- 
able, because the pseudogap appears to be caused by 
AF fluctuations! Third, in their recent paper Tallon and 
Loram argued, based on experimental results, that 
T* falls from a high value onto the Tc line instead of 
smoothly merging with Tc in the slightly overdoped re- 
gion. The above scenario for the pseudogap predicts the 
latter behavior of T* . In this paper we demonstrate that 
induced local spin-singlet amplitude due to short-range 
spin correlations causes a normal state pseudogap with 
d-wave symmetry even in the absence of pairing interac- 
tions. 

First of all we argue that there are two energy scales 
in the problem, because the pseudogap appears as a 
crossover phenomenon according to experiments. The 
low energy (or long distance) physics of AF and SC 
correlations is well captured by a static mean-field ap- 
proach, while the relatively high energy (or short dis- 
tance) physics of the pseudogap is invisible in such a 
study. Thus we resort to fluctuation theory in order to 
describe the dynamical nature of the pseudogap, and to 
determine T* and Apg. Note that the present formula- 
tion below is different from standard fluctuation theory 
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based on the (static) mean-field state. The mean-field re- 
sult of the i — J model will be used below solely to find the 
onset of leading correlations , and to compute mean- 
field AF and SC order parameters for the calculation of 
local spin and spin-singlet amplitudes. The mean-field 
t — J Hamiltonian reads 

Hmf ^Y.^{k)c\^^c^^^ 

k.{T 

— 2Jmy^(ct - cr 1^ — ct - c,- i ) 

k 

-Js^0d(fc)(ct^c^^ +cg^^c_,, p , (1) 
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where e{k) ~ — 2te(cosfca; 4- cosfcj,) — /i with x the hole 
density, m and s are mean-field AF and SC order pa- 
rameters determined from 
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FIG. 1. (a) Calculated mean-field phase diagram in doping 
{x = 1 — n) and temperature (T) plane in the t — J model for 
t/J = 3.0. T^''' and T^'' are mean-field AF and SC ordering 
temperatures, respectively. The filled diamonds are the pseu- 
dogap temperature determined from the single particle spec- 
tral function, (b) Interaction induced local spin-singlet (solid 
curve) and spin (dashed curve) amplitudes for t/J — 3.0 and 
T = 0.2J. 



where N is the total number of lattice sites. 4>d{k) = 
cos kx — cos ky and Q is the (commensurate) AF wave 
vector (tt, tt) in two dimensions. In this paper we restrict 
ourselves to a uniform solution which is just enough for 
our purpose. In a mean-field approximation, mean-field 
order already sets in when the correlation length reaches 
roughly one lattice spacing. This forces the above mean- 
field phase line (Fig. |^(a)) to be interpreted as the on- 
set of the corresponding short-range correlations. We 
identify Tj^^^ with another crossover temperature at 
which some magnetic experiments such as Knight shift 
show their maximum. For the parameter {t/J — 3.0) 
used in this paper, short-range spin correlations disap- 
pear at a; = — 0.19 — 0.20 at low temperature. 

We introduce AF and spin-singlet ||l^ correlation func- 
tions 

Xsp{i,T) = {Tra^{i,T)a^{0,0)) , 

Xpp(z,r) = (T,A<j(*,T)A^(0,0)), (3) 

where az{i) = (cl^c^j - c- ^Ci^i) = 2Sz{i) and Ad{i) = 
lJ2s9{S){ct+s,'\Ci^i - Ci+s,iCi,^) with g{S) = 1/2 for 6 = 
(±1,0), -1/2 for S = (0,±1), and otherwise. The AF 
and spin-singlet correlation functions are related to the 
local spin and spin-singlet amplitudes through the sum 
rules 
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where q = {q, ivm) and z/^ is bosonic Matsubara frequen- 
cies. T is absolute temperature. These sum rules can be 
easily obtained by taking r = 0^ limit and setting i to 
the origin in Eq. (|). In ter ms of renormalized vertices 
Usp and Vpp, [ p^ we approximate the AF and spin-singlet 
correlation functions 
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where the irreducible susceptibilities are defined as 
^EGO(fc-(z)GO(fc), 



Xphil) 



N ■ 



= ^ E(<^^(^) + '^<^^^~ k)fG\q k)G\k) . (6) 
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G^{k) is the noninteracting Green's function obtained 
from Eq. (Q) with J = 0. Now two unknown vertices, 
Usp and Vpp, are determined by the sum rules Eq. (|). 
|13| Hence, an increase in the local spin or spin-singlet 
amplitude evaluated in the interacting state over that in 
the noninteracting state leads to a nonvanishing positive 
value of Usp or Vpp , namely, the enhancement of the cor- 
responding correlation function. This (non-perturbative 
sum rule) approach has been shown to be quite reliable 
jlSf as long as short range correlations are concerned. In 
our calculations, the pseudogap appears when the spin- 
singlet correlation length reaches about 1 lattice con- 
stant. Since this method is expected to be accurate up 
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to the intermediate coupling regime, it forces the effec- 
tive bandwidth {W = 8tx) to be larger than the effective 
interaction strength (2J), which leads to a; > 0.08. In or- 
der to determine whether the pseudogap is caused by the 
spin-singlet or AF spin fluctuation channel, we separately 
consider the self-energy coming from each channel 

T 



J^spik) = UUsp- J2 Xsp{q)G\k - q) 



q 

T 
'N 
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{Mk) + k) fx„{q)G''{q - k) , (7) 



where U = 2J and V ^ J from Eq. (|l]). 

First let us begin by showing the interaction-induced 
local spin (dashed curve) and spin-singlet (solid curve) 
amplitudes (Fig.|l](b)) evaluated in the mean- field state 
of the t — J Hamiltonian in a region where s = (or 
T > T^'^). Since s = 0, the spin-singlet amplitude 
is entirely caused by short-range spin correlations in 
the absence of pairing interactions. Although in gen- 
eral a mean-field state is not accurate for strongly cor- 
related electron systems, certain local and short-range 
static quantities such as double occupancy and nearest 
neighbor correlations are reasonably well captured by the 
mean-field state particularly with AF order (See Ref. fl^ ] 
for more details). In fact the interaction-induced local 
spin and spin-singlet amplitudes (Eq. (^) arc determined 
most crucially by these quantities. The local spin ampli- 
tude starts to appear when short-range spin correlations 
begin to develop and keeps growing with decreasing dop- 
ing, as can be easily expected from the mean-field phase 
diagram (Fig.|l|(a)) itself. Quite unexpectedly, however, 
the local spin-singlet amplitude also increases with de- 
creasing doping despite the fact that the mean-field SC 
order s is absent. The increase of local spin-singlet am- 
plitude traces back to the growing short-range spin cor- 
relations with decreasing doping. This same feature was 
recently studied by the author in the context of the 
Hubbard model. 

In Fig. ^ we show the spectral functions at k = kp 
along (0, 0) — (0, tt) direction from spin-singlet (solid 
curve) and AF spin fluctuation (dashed curve) channels 
for two doing levels x = 0.15 and x — 0.08. For both 
densities the pseudogap appears first in the spin-singlet 
channel. This is verified even close to half-filling by 
considering the Hubbard model {U = 8t) [Q, which is 
not shown in this paper. The reason why the pseudogap 
is always caused by the spin-singlet channel is that as the 
local spin amplitude increases with decreasing doping, 
the local spin-singlet amplitude also increases at the same 
time. Hence the feature found away from half-filling per- 
sists down to half-filling. In our calculations the pseudo- 
gap appears when the characteristic low frequency scale 
of the spin-singlet correlations is smaller than tempera- 
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FIG. 2. Single particle spectral functions from spin-singlet 
(solid curve) and AF spin fluctuation (dashed curve) channels 
for (a) X = 0.15, T = 0.15J and for (b) x = 0.08, T = 0.41J. 



ture (renormalized classical regime). |13 The pseudogap 
due to AF spin fluctuations starts to appear for x < 0.10 
with r < T*. 

Based on the above results we show the pseudogap 
temperature T* (filled diamonds) as a function of dop- 
ing in Fig. |l|(a). T* falls from a high value onto the 
Tc (< T^"") fine instead of sharing a common line with 
Tc in overdoped region. When superconductivity is sup- 
pressed by setting s — 0, T* vanishes near Xc where 
short-range spin correlations disappear. It is not surpris- 
ing to find that T* closely follows T^^^ = T°, because in 
our study the pseudogap is caused by induced local spin- 
singlet amplitude due to short-range spin correlations, 
which is reasonably well captured by the mean-field state 
with AF order. All these features are at least qualita- 
tively consistent with the findings by Tallon and Loram. 
pO| Although quantitative agreement with the Hubbard 
and t ^ J models is achieved only for very strong cou- 
pling {U ^ < or J <C t), it is instructive to compare 
the calculated T* with the recent calculations obtained 
from the dynamical cluster approximation (DCA) for the 
Hubbard model by Jarrell et al. jl|] . Their T* near half- 
fiUing is about 0.09i for U = 8t, while ours is 0.15 - 0.16i 
for U = At^/J = 12t. The reasonable agreement with 
the more systematic approach is encouraging in light of 
the drastically simple approximation used in the paper, 
namely, replacing the strongly correlated hopping term 
of the t — J model by tx and using the mean-field state 
to compute local correlations. 

Figure ||(a) shows the pseudogap size Apg (filled di- 
amonds) which is defined as half of the pcak-to-peak 
distance in Fig. || by setting s = at T = 0. In the 
same figure the pseudogap energy extracted from vari- 
ous experiments by Tallon and Loram ||l^ is also shown 
as empty symbols for comparison . Apg vanishes near 
Xc, suggesting the presence of a quantum critical point 
at a critical doping. The agreement between our results 
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FIG. 3. (a) Calculated pseudogap size at T = denoted as 
filled diamonds for t/J — 3.0 and J — 125 meV. The open 
circles, diamonds, squares, up-triangles, and down-triangles 
are the pseudogap size extracted from susceptibility, heat ca- 
pacity, ARPES, NMR, and resistivity measurements, respec- 
tively, by Tallon and Loram. [|l^ (b) Total excitation gap 
Atg (open diamonds) and pseudogap Apg (filled diamonds) 
at T = 0. The inset shows their relative ratio. 



and experiments appears remarkable for such a simple 
approximation. The hnear vanishing of Apg near Xc is 
closely related to the corresponding behavior of the in- 
duced local spin-singlet amplitude. 

The total excitation gap (or ARPES leading edge gap 
or SC gap) in the SC state, Atg, can be also calculated. 

is always larger than Apg at T = due to 



Atg at T 

the additional contribution to the local spin-singlet am- 
plitude from s 7^ 0, which is shown in Fig. |(b). Apg, Atg 
and their relative ratio Apg/ Atg should be all monotoni- 
cally decreasing functions of doping, as shown in the inset 
of Fig. 0(b). Since the SC order parameter vanishes at 
Tc (at in this paper), the SC gap below Tc contin- 

uously evolves into the normal state pseudogap above Tc 
with the same momentum dependence and magnitude. 

In order to confirm the present scenario for the pseu- 
dogap, we propose an experiment in which long-range su- 
perconductivity is completely destroyed by a phase sen- 
sitive external perturbation such as a strong magnetic 
field. Or this can be also done in a vortex core with 
a relatively weak magnetic field. Then our scenario pre- 
dicts that the underlying ground state will manifest itself 
as an insulator with a pseudogap or spin gap. The nor- 
mal state pseudogap Apg should be observable down to 
T = for X < Xc vanishing near Xc, and its size is given 
by Fig. ^(a). The present scenario for the pseudogap pre- 
dicts that a normal state pseudogap is likely to appear 
when short-range spin correlations are well established 
and are not masked by long-range (AF or SC) order. 

The pseudogap appears here only as the suppression of 
low frequency spectral weight in certain physical quan- 



tities which are obtained through A(fc, w) or its convo- 
lution with a relevant vertex. It does not appear as a 
thermodynamic phase with broken symmetry. The pseu- 
dogap obtained in this paper is different from the spinon 
gap found in the previous slave boson [static) mean-field 
study of the t — J model. ||l8| In fact the pseudogap size 
and temperature obtained in the latter are the same as 
2Js and T^'^ calculated from Eq. |l| with m = 0. 

In the low energy sector short-range AF correlations 
compete with SC correlations to make or Tc go 

to zero near half-filling. At the same time, in the rela- 
tively high energy sector of order of J (or in the short 
distance scale) the same AF correlations induce the lo- 
cal spin-singlet amplitude, which is responsible for the 
normal state pseudogap in the present scenario. When 
the spin-singlet (or AF spin fluctuation) aspect is com- 
pletely neglected, the current approach reduces to the 
AF (or SC for T < T*^^) fluctuation scenario for the 
pseudogap. 0(2) SC fluctuations associated with super- 
fluid stiffness come into play below T*^^ |^ instead of 
T*, and diverge at Tc- The present results are robust to 
variations of t/J ~ 2.5 — 3.5 and small to moderate 
value of t' . 

In this paper we have considered only the local spin- 
singlet amplitude induced by short-range spin correla- 
tions and its consequences. The complete theory of long- 
range d-wave superconductivity is beyond the scope of 
the present formulation. How local spin-singlets (appear- 
ing at T^^^) acquire local SC phases (at Tc ^^) and even- 
tually establish their long-range phase coherence (at Tc) 
is a challenging problem to the theory of high tempera- 
ture superconductivity. 

In summary, we have shown that the local spin-singlet 
amplitude with d-wave symmetry, (|Ad(0)p), can be in- 
duced by short-range spin correlations even in the ab- 
sence of pairing interactions. Fluctuation theory has 
been formulated to make connection between T*, Apg 
and (|Ac((0)p). In the present scenario for the pseudo- 
gap, the normal state pseudogap is caused by the induced 
local spin-singlet amplitude due to short-range spin cor- 
relations, which compete in the low energy sector with 
SC correlations to make Tc go to zero near half-flUing. 
Since the SC order parameter vanishes at Tc (at T*^^ in 
this paper), the SC gap below Tc is smoothly connected 
to the normal state pseudogap above Tc with the same 
momentum dependence and magnitude. Calculated T* 
falls from a high value onto the Tc line and closely fol- 
lows T*^^. The calculated Apg is in good agreement with 
experimental results. We have proposed an experiment 
in which the present scenario can be critically tested. It 
would be interesting to see how robust are the features 
found in this paper, when the no-double-occupancy con- 
straint is strictly imposed on the t — J model and an 
inhomogeneous solution is used. 
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